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We present a semiclassical description of the level density of a two-dimensional circular quantum 
dot in a homogeneous magnetic field. We model the total potential (including electron-electron 
interaction) of the dot containing many electrons by a circular billiard, i.e., a hard- wall potential. 
Using the extended approach of the Gutzwiller theory developed by Creagh and Littlejohn, we 
derive an analytic semiclassical trace formula. For its numerical evaluation we use a generalization 
of the common Gaussian smoothing technique. In strong fields orbit bifurcations, boundary effects 
(grazing orbits) and diffractive effects (creeping orbits) come into play, and the comparison with the 
exact quantum-mechanical result shows major deviations. We show that the dominant corrections 
stem from grazing orbits, the other effects being much less important. We implement the boundary 
effects, replacing the Maslov index by a quantum-mechanical reflection phase, and obtain a good 
agreement between the semiclassical and the quantum result for all field strengths. With this 
description, we are able to explain the main features of the gross-shell structure in terms of just 
one or two classical periodic orbits. 
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1 INTRODUCTION 



1 Introduction 

The two-dimensional free-electron gas (2DEG) that 
occurs at the interface of suitably designed semicon- 
ductor heterojunctions (see, e.g. has attracted 
a lot of interest in the last years. This is partly due 
to the two dimensionality, which gives rise to new 
physical effects (such as the quantum Hall effect) 
and partly to the extremely high mobility of the 
electrons, which comes from the absence of (scatter- 
ing) donors or acceptors in the plane of the electron 
gas. The most attractive feature, however, is the 
great variability of these systems. With electron- 
beam lithography additional lateral constraints of 
the 2DEG down to structures of some 10 ran can 
be realized. This length is well below the typical 
phase coherence length and the electron mean free 
path (in GaAs both of them can be of the order 
of some jum), and can even be comparable to the 
Fermi wavelength of the electrons (typically 40 nm 
for GaAs), so that in such structures quantum con- 
finement effects play an important role. These sys- 
tems are therefore accessible on a quantum scale, 
opening up tremendous new possibilities in device 
design. 

Various approaches heave been used to model 
the 2DEG with and without additional lateral con- 
finement. Quantum-mechanical calculations tend 
to be rather involved and even for the simplest sys- 
tems numerically very demanding. Classical ap- 
proaches have the severe drawback that they ig- 
nore quantum interference effects, and are there- 
fore applicable only if the system dimensions are 
long compared to the mean free path and the phase 
coherence length. In the resulting gap of the theo- 
retical description, semiclassical approaches appear 
very promising. They approximate quantum me- 
chanics in such a way that the quantities involved 
can be interpreted classically, often in terms of the 
classical orbits in the system. They combine the 
advantages of the classical description, especially 
its limited numerical demands, with the ability to 
reproduce quantum-mechanical interference effects. 
This makes semiclassical methods a very attractive 
tool for mesoscopic systems, i.e., systems with di- 
mensions comparable to the phase coherence length 
and the mean free path. One of the most striking 
successes in recent years has been the explanation of 
conductance oscillations in superlattices, the Weiss 
oscillations Q. 

In this paper, we consider the level density of 
a 2DEG confined by external electric fields to a 



circular domain, with an additional homogeneous 
magnetic field perpendicular to the plane of the 
2DEG. When this quantum dot contains many elec- 
trons, the effective single-particle potential (i.e., the 
Kohn-Sham potential in the language of density 
functional theory, which contains the electron-elec- 
tron interaction in the local density approximation) 
is Wood-Saxon-like, with a flat region in the inte- 
rior and a rather steep surface. [] The level density 
is not too sensitive to details of the potential edge, 
so that a circular disk with infinite reflecting walls, 
i.e., a circular billiard, is a realistic model. 

The level density itself is hard to access exper- 
imentally, but it enters in many observable quan- 
tities. Persson et al. ||, for example, consider a 
quantum dot that is connected by two point con- 
tacts to the surrounding 2DEG. They propose an 
approximation in which the conductivity of this sys- 
tem in weak external magnetic fields is proportional 
to the level density of the dot at the Fermi en- 
ergy.0 Their measurements on a circular dot with 
about 1000 - 1500 electrons in a homogeneous mag- 
netic field show characteristic conductance oscilla- 
tions that could be well explained qualitatively in a 
perturbative approach by Reimann et al. J7|. They 
reproduce the oscillations in a simple and intuitive 
way by a few classical periodic orbits of the sys- 
tem and the flux enclosed by them. Because of its 
perturbative nature, this description only holds in 
weak fields. Another example where the level den- 
sity enters observable quantities is the magnetiza- 
tion IJ. 

For the interpretation of experiments on circular 
quantum dots, a semiclassical approximation of the 
level density with arbitrary field strength is desir- 
able. Such a description in terms of classical orbits 
is also of theoretical interest. In the absence of a 
magnetic field the classical orbits consist of straight 
paths bouncing at the boundary (see Fig. |2|). They 
have a one-dimensional degeneracy corresponding 
to the rotational symmetry of the system. In very 
strong fields, the confinement is negligible and the 
level density is dominated by the quantization of a 
free-electron gas, leading quantum mechanically to 
the Landau levels and described semiclassically by 
closed cyclotron orbits. These orbits do not touch 

1 This has been shown in self-consistent calculations for 
quantum dots || and is analogous to the situation in three- 
dimensional metal clusters Q. 

2 Recent exact quantum-mechanical calculations of the 
transport properties, however, show that this approximation 
is only valid for contacts at the opposite sides of the dot II . 



3 



the boundary and have a two-dimensional transla- 
tional degeneracy. A unified semiclassical descrip- 
tion thus has to include the transition between these 
two limiting cases, which includes changes of the 
topology and the degeneracy of the classical orbits. 
Treating these is of conceptual interest since both 
effects are well known to lead to divergences in semi- 
classical theories. 

We conclude this Introduction with a short out- 
line of the paper. As a reference, we first present in 
Sec. § the quantum-mechanical solution for the cir- 
cular billiard in homogeneous magnetic fields. Sec- 
tion H gives a short introduction to semiclassical 
methods, and in Sec. || we derive a semiclassical 
trace formula of the disk. We then compare its 
results to the quantum-mechanical ones for vari- 
ous field strengths. The agreement is good for very 
weak and for very strong fields, but the semiclas- 
sical approach rather appears to fail in the inter- 
mediate regime. The deviations are due to bifurca- 
tions of classical orbits, to diffraction effects, and to 
boundary effects. The latter give the largest contri- 
butions, and in Sec. || we develop a simple approxi- 
mation to include these effects in the trace formula. 
This corrected trace formula gives satisfactory re- 
sults for all magnetic field strengths, and we give 
an intuitive interpretation of the the level density 
in the various -B-field regimes. The paper closes 
with a summary of the results and an outlook to 
further investigations. 



where 

i*i ^-a n vA + \i\;^] =0 . (2) 

Here n > denotes the radial and I the angular- 
momentum quantum number. The zeros of \F\ 
were determined numerically — which is concep- 
tually easy but requires a lot of numerical work. 
Figure [I] shows the well-known dependence of the 
eigenvalues E n \ on B. One clearly sees how with in- 
creasing magnetic field the different states condense 
into the Landau levels (dashed lines). 
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Figure 1: The quantum- mechanical eigenenergies of 
the circular billiard in dependence of the magnetic field. 
The dashed lines correspond to the four lowest Landau 
levels. 



3 Semiclassical methods 



2 The quantum-mechanical solu- 
tion 

In the following, we will use normalized energies 
E in units of h 2 /2mR 2 and normalized magnetic 
fields B in units of h/eR 2 , where R is the disk ra- 
dius. In these units, we have \f~E~ = kR and with 
the classical cyclotron radius R c = Hk/eB, we get 
R c /R = kR/B. 

The exact quantum-mechanical solution for the 
circular billiard in homogeneous magnetic fields was 



presented by Geerinckx [1C] and, using a different 



approach, by Klama and RoBler [11]. The eigenen- 
ergies are given by the zeros of the confluent hyper- 
geometric function \F\ as 
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For semiclassical approximations, a broad variety of 
methods is at hand. Some of them approximate di- 
rectly the quantum-mechanical eigenvalues |[^, [b| . 
whereas others describe the level density g(E). For 
this purpose, g(E) is usually split up into a smooth 
part g(E), the (extended) Thomas-Fermi level den- 
sity, and an oscillating part 5g: 



g(E) = g(E) + Sg(E) . 



(3) 



The latter can be expressed in terms of the periodic 
orbits of the corresponding classical system, which 
is therefore called periodic orbit theory (POT). Such 
relations have been established by various approaches 
[14, 15, [l6|, [I?], usually resulting in so-called 



trace formulae of the form 
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4 TRACE FORMULA FOR THE CIRCULAR BILLIARD 



Here T labels all classical periodic orbits of the sys- 
tem. Each orbit contributes to the level density via 
an oscillating term that depends on the classical ac- 
tion Sr along the orbit and on the Maslov index or 
depending on the orbit's topology. The amplitude 
A-p is a slowly varying function of energy, deter- 
mined by classical properties of the orbit such as 
its degeneracy and its stability^] All semiclassical 
approaches have their individual merits and draw- 
backs |l{J, and it is interesting to note that for the 
simple integrable case of the circular billiard with- 
out magnetic field, all applicable methods result in 
the same trace formula [20, |2l|, ^2|, which further- 



more reproduces exactly |22] the EBKQ spectrum. 
Applying a homogeneous magnetic field, the sys- 
tem remains integrable, with the energy and the z 
component of the conjugate angular momentum as 
constants of the motion. In weak fields, this sys- 
tem was treated using a perturbative approach by 
Bogachek and Gogadze [21], Ullmo et al. [j|], and 
by Reimann et al. Q. We have chosen the trace 



formula of Creagh and Littlejohn jl8| starting 
point for the description in arbitrarily strong fields. 



4.1 Classification of the periodic orbits 

The complete classification of the periodic orbits of 
our system is straightforward. Let us first consider 
the case without magnetic field. In a circular bil- 
liard, the periodic orbits (PO) are identical to those 
in a three-dimensional spherical cavity, whose com- 
plete classification has been given by Balian and 
Bloch All orbits have a one-dimensional de- 

generacy corresponding to the rotational symmetry 
of the system. Each family of degenerate orbits 
with a given action (or length) can be represented 
by a regular polygon. The first few polygons are 
shown in Fig. 0. These orbit families are classi- 
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Trace formula for the circular 
billiard 



The trace formula of Creagh and Littlejohn [|18[ is 
well suited for the semiclassical description of the 
circular billiard, as it can deal with continuous sym- 
metries. The main idea of their approach is the 
separation into a symmetry-free system treated by 
usual semiclassical techniques and the symmetry, 
which is used to integrate over the orbit families. 
The structure of the trace formula (Q) remains es- 
sentially unchanged by this procedure, but the def- 
inition of A is different, reflecting the different clas- 
sical structure of the dynamics. For the details we 



refer to the original publication [18]. In order to 
calculate the level density with this trace formula, 
we have to classify the periodic orbits and calcu- 
late their actions, amplitudes, and Maslov indices. 
These steps are presented in the following subsec- 
tions. 



3 Especially, A depends on h only by a factor h~ k ^ 2 , where 
k is the degree of degeneracy of the orbit family. 

4 EBK stands for the semiclassical approximation devel- 
oped by Einstein, Brillouin and Keller 113] 



Figure 2: The classical periodical orbits of the circular 
billiard in the absence of a magnetic field are the regu- 
lar polygons. They can be classified with (v,w), where v 
is the number of corners and w indicates how often the 
trajectory winds round the center of the disk. 



fied [15| by /3 = (v , w), where v denotes the number 
of corners (vertices), and w is the winding number, 
i.e., it counts how often an orbit winds around the 
center. With v > 2w > 2, (3 = (v, w) uniquely 
describes all families of POs of the system in the 
absence of a magnetic field. Because of the time- 
reversal symmetry, all orbits except the diameter 
(v = 2w) have an additional discrete twofold de- 
generacy, which has to be accounted for in the trace 
formula. 

Switching on the magnetic field causes the clas- 
sical trajectories to bend, with the direction of the 
curvature depending on the direction of motion with 
respect to the magnetic field. This entails a break- 
ing of time-reversal symmetry. For weak fields, the 
orbits can still be classified by (3 if an additional 
index (±) is introduced. This situation is shown 
in the upper row of diagrams in Fig. [3| for the orbit 
= (4, 1). Up to a field strength where the classical 
cyclotron radius R c equals the disk radius R, hence- 
forth referred to as the weak-field regime, the orbits 
do not change their topology and the classification 



4.2 The bouncing orbits 
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f3 ± holds. For the strong-field regime with B > kR, 4.2 The bouncing orbits 



the structure of the POs is different. This situation 
is shown in the second row of diagrams in Fig. ||. 
The (3~ orbits change their shape continuously over 
the point R c = R, but the f3 + orbits change their 
topology abruptly. However, since there is a one- 
to-one correspondence between orbits for R C >R and 
for R C <R, still gives a complete classification of 
all bouncing orbits, i.e., of orbits that are reflected 
at the boundary. For R c < R, there are additional 
cyclotron orbits that do not touch the boundary at 
all. They have to be included separately in the 
sum over all orbits in the trace formula. At field 




(4,1) 
Rp < R : 



(4,1) 




(4,1)" 



(4,ir 



cyclotron orbits 



Figure 3: A magnetic field breaks the time-reversal 
symmetry, so that the orbits are no longer independent of 
the direction of motion. Introducing an additional index 
±, the orbits can be classified by (v, w) , both in weak 
(R c > R) and in strong (R c < R) fields. For strong 
fields, there occurs an additional family of orbits that do 
not touch the boundary, the cyclotron orbits. 

strengths where R c < Rsm{irw /v ), the (vjw)^ or- 
bits no longer exist (see Fig. |j). They vanish pair- 
wise, which is the simplest case of an orbit bifur- 
cation. This imposes an additional restriction on 
the sum over (v,w). Including this, we now have a 
complete classification of all periodic orbits in the 
circular billiard at arbitrary field strengths. 



(v,w) + 






(v,w) 



R c > R sin(7t w/v) R c = R sin(7t w/v) R c < R sin(n-w/v) 

Figure 4: At afield strength where R c = i?sin(7r w/v), 
the orbits (vjw)^ vanish pairwise. 



The action of a closed orbit in a magnetic field can 
be written as the sum of the kinetic part and the 
magnetic flux enclosed by the orbit 



S, 
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pdq = hkLp — eBFp 



(5) 



The geometrical lengths Lp and the enclosed areas 
Fp of the periodic orbits discussed above (correctly 
counting those areas that are enclosed several times, 
cf. Fig. [5]) can be calculated by elementary geome- 
try. In terms of the geometrical quantities R c ,R,^f 
and O, explained in Fig. O, we obtain 






(3,l,n) n 



(5,2,n) n 



(7,3,n) n 



Figure 5: Calculating the magnetic flux enclosed by an 
orbit, the multiply enclosed areas (darker gray) have to 
be correctly accounted for. 
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According to the trace formula [18], the orbit am- 
plitudes are composed of an integral over the sym- 
metry group, which for the rotational U(l) sym- 
metry of the disk just gives 2tt/v, of the period of 
the orbit L/fik, and of the Jacobian resulting from 
the symmetry reduction dL/d^>, where \& = — 2n0. 
All these quantities can be calculated analytically, 
resulting in 



4> 
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1 




E y/Rk^ ^ R^ sR 

-7 for (f3 + ,R c < R) 
otherwise 



(7) 



where c, d, and s are the geometrical lengths sketched 
in Fig. ^. The connection of these geometrical quan- 
tities to the classification parameter [3^ and the cy- 
clotron radius R c is given in Appendix |A|. 
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4 TRACE FORMULA FOR THE CIRCULAR BILLIARD 




Figure 6: The actions and amplitudes of the classical 
periodical orbits can be expressed purely in terms of the 
geometrical quantities shown here. 

4.3 Cyclotron orbits 

For magnetic fields stronger than B = kR, the clas- 
sical cyclotron radius R c is smaller than the disk 
radius R. This gives rise to a new class of periodic 
orbits, the cyclotron orbits, which do not touch the 
boundary at all (see Fig. [|). They form transla- 
tionally degenerate families, whereas the bouncing 
orbits (v, w)^ considered above are degenerate with 
respect to rotations. For the translational case, 
the symmetry reduction can be performed directly, 
without need of the general procedure of Creagh 
and Littlejohn. We transform the phase-space co- 
ordinates according to 



1 



n x := 



eB 



yrr\ \ Px + ~Y y ' ' n * := % + \/l eB l x 



i 



Py ~ -ir x ) ; u y '■= 



eB\ y . 

(8) 




Figure 7: The motion of a charged particle in a ho- 
mogeneous magnetic field can be expressed in the coordi- 
nates of the relative motion (x,y) — |eS| _1 / 2 (— ir y , ir x ) 
and the coordinates of the center of gyration (X, Y) = 
\eB\~ 1 / 2 (Jl x ,—I[y). The Hamiltonian is independent of 
(H x ,Hy); all orbits with the center (X, Y) in the gray 
shaded area are degenerate. 

accessible for (11^,11^), which can be directly read 
off Fig. ^ (shaded area). We therefore get for the 
degeneracy 
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(10) 



Now, the Hamiltonian Eq. (||) is identical to that 
of a one-dimensional harmonic oscillator. Using its 
analytically known trace formula,^] the contribution 
of the cyclotron orbits to the oscillating part of the 
level density is given by 



6g c 
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2 oo 



cos(nkirR c 



nix 
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(11) 

Here n is the winding number around the center of 
gyration. Note that the frequency is again deter- 
mined by the classical action along the orbit, which 
in this case is 



Apart from a factor ^/|eS|, (ir x ,ir y ) are the coor- 
dinates of the motion relative to the center of gy- 
ration (IL x ,H y ), as illustrated in Fig. |7[ In these 
coordinates the Hamiltonian reads 



S = n ■ hk ■ ttR c 



(12) 



Note that here exactly half of the kinetic contribu- 
tion to the action is canceled by the flux term. 



H 



eB 
2m 



( 9 ) 4.4 Additional phases 



x and Hy are 

[n s ,n y ] = 



As expected, H does not depend on the coordi- 
nates of the center of gyration. H 
canonically conjugate variables, since 
ih. Because the relative and the center-of-gyration 
coordinates commute, i.e., [n^,^] = [n^, TT y ] = 
[Hy,ir x ] = [YL y ,ir y ] = 0, the degeneracy of a cy- 
clotron orbit is simply the phase-space volume V 



For a discussion of the additional phases in the trace 
formula (f|) we refer to Sec. 5.1. There we find that 
the Maslov index for bouncing orbits is a = 3v, 
and for the cyclotron orbits it is a = 2. According 



to [18, 24 1, we have an additional phase of 5 tt/2 



5 The harmonic oscillator is one of the few cases that can 
be treated exactly within standard POT |23| . 



4.6 Results in the weak-field regime 



7 



stemming from the symmetry reduction. 5 is found 
to be 



for (P+,Rc < R) 

1 otherwise. 



(13) 



We have now analytic formulas for all quantities 
of the trace formula. The numerical evaluation of 
this semiclassical level density will be performed in 
Sees. O and|4~7. 



4.5 The shell structure 

In an experiment, the observed levels are always 
broadened due to temperature, life-time, or impu- 
rity effects. In most systems the levels are not ap- 
proximately equally spaced, but occur in bunches, 
the so-called shells, which are separated by rela- 
tively wide energy gaps. Smoothing the level den- 
sity over a width larger than the typical level spac- 
ing, but smaller than the distance of the bunches, 
reveals the (gross-) shell structure of the system. 
It contains in many cases the dominating quantum 
effects. This folding procedure can easily be im- 
plemented in the semiclassical trace formula. For 
pure billiard systems, i.e., systems where S = kL 
with L independent of k, a Gaussian folding of the 
level density is equivalent to the multiplication of 
the orbit amplitudes in the trace formula with a 
Gaussian with reciprocal width. In Appendix ^ we 
give a more general form of this relation, which is 
not restricted to pure billiard systems and to Gaus- 
sian smoothing, but can deal with general systems 
and arbitrary smoothing functions. We will in the 
following use this generalized approach, as in finite 
magnetic fields we no longer have a pure billiard 
system. As this is a more technical point, we leave 
the discussion for Appendix [l|. There we also give 
detailed information about the numerical evalua- 
tion scheme and the smoothing function used. The 
latter is in the following characterized by a parame- 
ter 7, which corresponds to the variance of a Gaus- 
sian exp [-1/2 (feR/7) ] wi t h the same half- width. 

The additional factor in the amplitudes stem- 
ming from the smoothing strongly suppresses the 
longer periodic orbits^] so that usually only a few 
of them (2 - 10) contribute to the gross-shell struc- 
ture. This makes the POT a very convenient tool 
for the calculation of this quantity. The quantum- 
mechanical approach is in some sense complemen- 
tary to the semiclassical one. It first gives the single 

6 This holds only for billiard systems. The suitable gener- 
alization is again given in Appendix [b| 



eigenvalues, of which many have to be known to cal- 
culate the shell structure. On the other hand, a full 
semiclassical quantization, i.e., resolving the level 
density down to the single eigenenergies, involves 
in general an exponentially increasing number of 
orbits and thus is a very demanding task. Here we 
are mainly interested in the semiclassical calcula- 
tion of the gross-shell structure, for which only a 
few of the shortest and most degenerate periodic 
orbits are required. We will, nevertheless, also try 
to go for a full quantization — mainly to verify the 
quality of our semiclassical approximation. 

4.6 Results in the weak-field regime 

In the previous sections we have derived an an- 
alytical trace formula for the circular billiard in 
homogeneous magnetic fields of arbitrary strength. 
Presently we shall discuss the resulting level densi- 
ties as a function of energy and magnetic field. Let 
us start with weak fields {R c > R), for which the 
topology of the classical periodic orbits is the same 
as in the absence of a magnetic field (see Sect. 4.1), 
so that we expect the semiclassical approach to be 
of the same quality as for zero field. The high-field 
regime, where we expect new effects to arise, will 
be the topic of the next section. 

The case of the circular billiard in small homo- 
geneous magnetic fields has already been treated 
by Bogachek and Gogadze |2l|] and by Reimann et 
al. [0] using a perturbative approach for weak fields. 
Replacing the amplitudes of Eq. ^ by their asymp- 
totic values for B —* and expanding the actions 
of Eq. (||) up to first order in B reproduces, indeed, 
their results. 

In Fig. |8] the semiclassical level density obtained 
with 7 = 0.35 (solid line) is plotted against the 
equivalently smoothed quantum result (dashed) for 
various values of B. The agreement is almost per- 
fect — just as it is in the zero-field case, which has 



been extensively discussed by Reimann et al. [25]. 
Note that the calculation requires over 850 numeri- 
cally determined eigenvalues for the quantum-mechani- 
cal calculation (which then have to be smoothed), 
whereas the semiclassical calculation is analytical 
and just requires the most important orbits (the di- 
ameter and the two triangle, square, pentagon and 
hexagon orbits). 

Since we have a classification of all periodic or- 
bits and analytic expressions for their actions and 
amplitudes, we can attempt a full semiclassical quan- 
tization by summing up sufficiently many of them. 
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Figure 8: The semiclassical level density of the disk 
billiard (solid line) matches perfectly the equivalently 
smoothed quantum-mechanical result (dashed, and well 
hidden under the solid line). The smoothing width is 
7 = 0.35. 

The result is shown in Fig. || where we display the 
total level density g = g + 6g, averaged over a width 
7 = 0.025, which is smaller than the typical level 
spacing. As Tanaka has shown in Q, the smooth 
part of the level density of the circular billiard does 
not depend on the magnetic field to leading order 
in h. We use the Thomas-Fermi level density for 
zero field, which is identical to the familiar Weyl 



expansion [26] 
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(14) 



Both the semiclassical level density (solid) and the 
corresponding quantum-mechanical one (dashed) in 
Fig. H exhibit clearly separated peaks whose heights 
give the degeneracies of the individual levels. The 
two lines can hardly be distinguished, thus the semi- 
classical approach gives almost perfect results even 
in this extreme case of full quantization. 

4.7 Results in the strong-field regime 

Figure is the strong-field equivalent of Fig. ||. 
It displays again the semiclassical (solid) and the 
quantum- mechanical (dashed) level densities, ob- 
tained with an equivalent averaging width 7 = 0.35. 
The agreement for small fields (R c < R) is good, as 



already shown in Sec. 4.6. For stronger fields, the 



positions of the Landau levels (gray lines in Fig. 10) 
are well reproduced, but their degeneracies are over- 
estimated in the semiclassical approximation. 

Figure 11 corresponds to Fig. |9| and displays 
the full quantization of the system. The semi- 
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Figure 9: The semiclassical level density (solid) for 
a smoothing width 7 = 0.025, which is small enough 
to resolve the single eigenenergies ("full quantization"). 
One can hardly distinguish this from the equivalently 
smoothed quantum-mechanical result ( dashed line under- 
neath). The positions of the quantum-mechanical eigen- 
values in dependence of B are indicated by the gray lines. 
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Figure 10: The semiclassical level density (7 = 0.35,) 
of the disk billiard (solid) compared to the equivalently 
smoothed quantum-mechanical result (dashed). The gray 
lines and the arrows indicate the positions of the first 
four Landau levels. For strong fields (R c < R), the 
agreement between the semiclassical and the quantum- 
mechanical results is not satisfactory. 

classical approach is seen to fail for stronger fields. 
As already mentioned, this is due to the neglect of 
various effects. First, there are orbit bifurcations 
where classical orbits vanish pairwise with increas- 
ing magnetic field (see Figs, f|, [L7] and |l8|). The 
change of the topology of the f3 + orbits and the oc- 
currence of cyclotron orbits are also bifurcation ef- 
fects. Those are known to lead to divergences in the 
trace formula. Second, we have neglected boundary 
effects from grazing orbits and diffraction effects, 
which could be implemented in the trace formula 
by considering creeping orbits. A closer look at 



Fig. 11 gives some hints as to which of these effects 
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Figure 11: The semiclassical level density resolved up 
to the single eigenenergies (solid) compared to the equiv- 
alently smoothed quantum-mechanical levels (dashed). 
The gray lines indicate the positions of the eigenener- 
gies as functions of B. For strong fields (R c < R) the 
agreement between semiclassical and quantum-mechani- 
cal calculation is bad. Note especially the large discrep- 
ancies in the degeneracies of the Landau levels, and the 
completely missing states slightly above the Landau levels 
(see insets). 

dominate. The two most striking observations are 
as follows: (1) The trace formula reproduces well 
the positions of the Landau levels,[| but it overes- 
timates the degeneracies of these states. The er- 
ror becomes smaller with increasing field strength 
(see insets in Fig. 0). (2) The levels that have 
energies slightly above the Landau levels are com- 
pletely missed by the semiclassical approach. These 
two observations suggest that it is a boundary ef- 
fect that causes the discrepancies. A simple hand- 
waving argument might be useful to illustrate the 
effect. Quantum mechanically, a particle moving 
on a cyclotron orbit will feel the boundary even if 
classically not touching it. Particles on cyclotron 
orbits close to the boundary thus feel an additional 
confinement. This restriction to a smaller volume 
will lead to a higher energy. In this picture, not 
all the cyclotron orbits are degenerate. The orbits 
close to the boundary will no longer have the energy 
of the Landau level, but a slightly higher one. This 
is exactly what would correct the observed defects 
of the semiclassical approximation. In the next sec- 
tion we will present a simple way to incorporate this 
boundary effect in the trace formula. 



7 This is no surprise since the Landau levels are due to the 
free cyclotron motion of the electrons, which is equivalent to 
that of a ID harmonic oscillator. The latter is known to be 
exact in the semiclassical approximation. 



The observations of Sec. 47 suggest that boundary 
effects are responsible for the failure of the semiclas- 
sical approximation in strong fields. The only place 
where boundary properties enter the standard trace 
formula is the Maslov index. We therefore propose 
here to replace the Maslov index by a more sophis- 
ticated quantity, which includes some quantum ef- 
fects. Before doing this, let us give a brief summary 
of the origin of the Maslov index. 

5.1 The Maslov index 

The origin of the Maslov index can most easily 
be understood in the one-dimensional case. Semi- 
classically, one approximates the wave functions to 
lowest order by plane waves with the local wave 
number k(x) = y / 2m[E — V(x)]. This approxima- 
tion obviously breaks down at the classical turn- 
ing points where E = V(x) and the wavelength 
diverges. Expanding the wave function around the 
classical turning points and matching the solutions 
to the plane-wave solutions far from the turning 
points leads to additional phases in the semiclas- 
sical quantization [12|. In the limit Ti — > these 
are independent of the detailed shape of the poten- 
tial. Each reflection at a softf^] turning point gives 
a phase of — vr/2, whereas each reflection at an in- 
finitely steep wall gives a phase of —it. Writing this 
phase as —a tt/2, one usually calls a the Maslov 
index. 

In the case of the circular disk, the Maslov in- 
dex can be obtained simply by counting the classi- 
cal turning points of the one-dimensional effective 
potential in the radial variable r. For skipping or- 
bits, the Maslov index per bounce is 3, including 
one soft reflection at the centrifugal barrier and one 
hard- wall reflection. For the cyclotron orbits, the 
effective potential is a one-dimensional harmonic os- 
cillator (see Sec. 4.3) with two soft turning points, 
and thus their Maslov index per period is 2. 

5.2 Reflection phases 

For finite h the additional phase stemming from 
classical turning points will depend on the shape of 
the potential. Let us consider a cyclotron orbit at 
a distance x\y to the billiard boundary. Neglecting 



"Soft" here means that the slopes of the potential at the 
classical turning points are finite. 
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the curvature of the boundary (which corresponds 
to the strong- field limit), we can reduce the mo- 
tion in the presence of the wall to an effective ID 
motion just as in the unbounded case presented in 
Sec. O. This is shown in Fig. The upper row 




Figure 12: The planar cyclotron orbit is equivalent 
to the motion in a one- dimensional harmonic oscillator 

(a) . Neglecting its curvature, the billiard boundary can 
be implemented in the effective one- dimensional motion 

(b) -(d). 

of diagrams shows the 2D motion, the lower row 
gives the reduction to the one-dimensional motion 
in an effective potential. Figure |12| (a) shows the 
unbounded case, in (b) the orbit is near the bound- 
ary, and (c,d) illustrate skipping orbits. 



hara and Ebina [27| who used local expansions in 
terms of Airy functions. We use a different ap- 
proach and integrate the quantum-mechanical prob- 
lem numerically. From the solutions we calculate 
the reflection phases (pn, which is displayed in Fig. 13 
As expected, they show a smooth transition from 




A particle in the potential sketched in Fig. |12[ b 
is classically not influenced by the additional wall, 
since it will never touch it. Quantum mechani- 
cally, however, the wave function enters the classi- 
cally forbidden region and thus feels the boundary 
even for xw > Rc- This leads to a smooth tran- 
sition of the quantum-mechanical reflection phase 
(fR over the point xw = R c , whereas the semiclas- 
sical Maslov phase is discontinuous at this point; 
as we have just seen in Sec. 5A above, it is — tt 
for xw > Rc and — 3/2ir for xw < R c . Our way 
to implement these quantum effects at the bound- 
ary in the semiclassical trace formula is therefore to 
replace the Maslov index by the quantum-mechani- 
cal reflection phase (fiR of the corresponding one- 
dimensional motion. This smooth version of the 
Maslov phase will also remove the former clear sep- 
aration between cyclotron orbits and skipping or- 
bits. These two limiting cases are now continuously 
linked, with ipR ranging between — tt and — 3/2n. 
We will refer to the orbits in the transition region, 
which are close to the boundary within h, as to the 
grazing orbits. 

The calculation of the reflection phases is in this 
approximation reduced to the problem of the one- 
dimensional harmonic oscillator in an additional square- 
well potential. This system was approached by Isi- 



Figure 13: The reflection phase ipn in dependence of 
the distance of the center of gyration from the bound- 
ary xw- The transition from xw < 1 to xw > 1 is 
continuous and gets sharper for increasing (kR) 2 / B . In 
the limit (kR) 2 / B —>■ oo, which corresponds to the semi- 
classical limit h — ► 0, the Maslov phase (thick line) is 
recovered. 

—tt at xw S> R c to — 3/27T at xw *C R c - The 
transition gets sharper if (kR) 2 /B increases. For 
(kR) 2 /B — > oo, which corresponds to the semiclas- 
sical limit h — > 0, the standard Maslov phase (thick 
line) is reproduced. Quantum corrections are seen 
to have the greatest influence on the grazing orbits 
(xw ~ Rc) and on orbits with xw~ — R c - These 
are known as the whispering gallery orbits, as they 
move in a narrow region along the boundary. 



5.3 Comparison to the 

quantum-mechanical result 

Figures |l4] and ^ show the coarse-grained level 
density and the full quantization of the spectrum, 
respectively, both calculated with the reflection pha- 
ses of Sect. |5,2|, A comparison with the corre- 



sponding diagrams in Figs. 1C and 11, which display 
the result obtained with the standard Maslov in- 
dices, immediately shows that the situation is dras- 
tically improved when using reflection phases. The 
coarse-grained level density now is good at all mag- 
netic field strengths. The full quantization displayed 



in Fig. 15 is not perfect, but the most striking error 



in standard POT, giving the wrong degeneracies of 
the Landau levels, is now corrected. This is dis- 
played in detail in the insets 1 and 2 in Figs. 11 



5.4 Semiclassical interpretation of the 



shell structure 
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Figure 14: 27ie semiclassical coarse-grained level den- 
sity of the disk calculated with reflection phases (solid) 
compared to the equivalently smoothed quantum-mechani- 
cal result (dashed). The gray lines and the arrows in- 
dicate the positions of the lowest Landau levels. The 
agreement is considerably better than with the use of the 
Maslov indices as displayed in Fig. uh. 
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Figure 15: The semiclassical level density with cor- 
rected reflection phases resolved up to the single eigenen- 
ergies of the billiard (solid) compared to the equivalently 
smoothed quantum-mechanical result (dashed). The gray 
lines indicate the positions of the eigenenergies in depen- 
dence of B. The agreement is much better than in the 
case of the Maslov indices in Fig. [71|. The degeneracies 
of the Landau levels are correctly reproduced, only the 
levels that are close to condensing on the Landau levels 
show deviations (insets). 



we have neglected here. The resulting error can be 
understood as follows: a generic two-dimensional 
system has two quantum numbers, thus requiring 
two semiclassical quantization condition. The free 
2D electron gas in a homogeneous magnetic field 
has an additional dynamical symmetry and only 
one quantum number (labeling the Landau level). 
This additional symmetry is broken by the pres- 
ence of a curved boundary.^ This implies that the 
semiclassical description in terms of cyclotron orbits 
near the boundary misses one quantization condi- 
tion, which is hidden in the broken dynamical sym- 
metry. Therefore these orbits give rise to a continu- 
ous semiclassical (sub) spectrum. Bouncing orbits, 
however, have the correct symmetry and lead to a 
discrete subspectrum. This transition can be seen 
in Fig. [D]. On the low-energy side of inset 3, the 
semiclassical level density shows a continuous spec- 
trum stemming from the grazing orbits, whereas 
the quantum-mechanical result gives quantized lev- 
els. This error affects mainly the fully quantized 
spectrum; the influence on the gross-shell structure 
is negligible. 

Figure [0] shows once again the semiclassical 
level density calculated with reflection phases, now 
in the whole range from zero field to full Landau 
quantization (solid). The comparison with the ex- 
act quantum result (dashed) shows that the semi- 
classical approximation is in fact valid for arbitrar- 
ily strong fields. Small deviations occur only at the 
bifurcation points of the dominating orbits. As al- 
ready mentioned, we did not include the effects of 
the bifurcations in our calculation. The resulting 
errors are much smaller than the effect of the re- 
flection phase, and they are seen to be more impor- 
tant for the gross-shell structure than for the full 
quantization.^ 

5.4 Semiclassical interpretation of the 
shell structure 



and 15, respectively. The single states between the 
Landau levels, however, are still not reproduced 
correctly. This is due to our simple approximation, 
which only includes boundary effects via the reflec- 
tion phase. The classical orbits are not changed, 
so that in our approximation the center of gyra- 
tion of the cyclotron orbits (x\y > Rc) is fixed, 
whereas for bouncing orbits (xw < Rc) it moves 
around the disk. Modeling the expected smooth 
transition from x\y > R c to xw < Rc semiclassi- 
cally would require including diffractive orbits that 



In Sec. 5.3 we have shown that the semiclassical 



approximation for the level density is valid for arbi- 
trarily strong fields. It reproduces the exact quan- 



A straight boundary does not break the symmetry. This 
is the reason why in this case it is possible to reduce the 
system to one dimension, which we have exploited in Sec. [l^ 
for the calculation of the reflection phase. 

10 This implies that even though the amplitudes are diverg- 
ing, the trace formula can still be used. Note, however, that 
near the bifurcation points the numerical evaluation of the 
trace formula has to be performed with special care, as de- 
scribed in Appendix H. 
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tum-mechanical result with a remarkably reduced 
numerical effort. For the quantum-mechanical cal- 



culation shown in Fig. 16 about 2500 eigenvalues 
had to be calculated and numerically smoothed for 
each value of B, whereas the semiclassical result is 
obtained summing the contributions of just 20 or- 
bits. The most attractive feature of the semiclas- 
sical approximation, however, is the simple, intu- 
itive picture it gives. Let us now exploit this to ex- 
plain the behavior of the shell structure of the disk 
billiard in terms of classical quantities. According 
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Figure 16: The semiclassical coarse-grained level den- 
sity of the disk billiard with corrected reflection phases 
(black) compared to the equivalently smoothed quantum- 
mechanical result (gray). The agreement is acceptable 
in the whole range of energies, disk radii, and magnetic 
fields. The vertical lines indicate the bifurcation points 
of the most important orbits. The shaded regions are en- 
larged in the figures below. There the thick lines show the 
interpretation of the level density as given in the text. 

to the trace formula Eq. (|j), each periodic orbit [3 
contributes an oscillating term to 5g. Its frequency 
is determined by the classical action Sp along this 
path, which can be locally approximated by 



Sp(k) = Sp(k ) + KGp(k) (k-k ) 



(15) 



with the quasiperiod TiG. As shown in Appendix |B|, 
for systems with constant absolute velocity along 
an orbit G is the geometrical orbit length. The 
amplitudes of the oscillating terms are ApF(Gp), 
where F is the window function that depends on 
the desired smoothing of the level density (see Ap- 
pendix [b]). Before we interpret the contributions of 
the various orbits to 5g, let us discuss the behavior 
of Gp and Ap. 

Figure |l7| shows the dependence^ of G on the ratio 
R c /R. Note that for R c > R (see right diagram 
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Figure 17: The quasiperiods G of the most important 
orbits in dependence of R c /R. For R c > R, G is inde- 
pendent of the index ±. The orbit bifurcation points in 
strong fields ( vertical lines ) can be clearly seen. 



of Fig. 17) G is independent of the direction of mo- 
tion ±, even if the classical action depends on it. At 
R = R c all orbits are creeping along the boundary, 
forming collectively the whispering-gallery mode. In 
strong fields (R c < R, left diagram) G is different 
for the "+" and the "-" orbits. Only at the bifur- 
cation points, where the two orbits coincide, they 
have identical G. For strong fields, the value of G 
at the bifurcation points converges to w ■ tt 2 R. 



In Fig. 18 the amplitudes of the orbits relative 
to the B = values, 



4° 





sin 



3/2 Q 



(16) 



are plotted versus the ratio kR/B = R c /R. The 
amplitudes of the "-" orbit is always larger than 
that of the corresponding "+" orbit. At R c = 
R, where the "+" orbits change the topology (see 
Fig. H), their amplitudes are zero, so that these 
bifurcations do not lead to artifacts in the level 
density. In stronger fields, the amplitudes diverge 
at the bifurcation points, indicating that the semi- 
classical approximation breaks down at these points 
(more exactly, one of the saddle-point approxima- 
tions in the derivation of the trace formula becomes 
invalid). A rigorous treatment of these bifurca- 
tions has been presented in a very general form 



For R c > R even 10 orbits are sufficient. 



2 The explicit formula for G is given in Appendix [b| 



5.4 Semiclassical interpretation of the 



shell structure 
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Figure 18: The amplitudes of the dominating orbits 
[3 = (i>, 1) with v — 2, . . . , 5 relative to their B = Q 
value. (The amplitude of the cyclotron orbit is in arbi- 
trary units.) At the bifurcation points R c = sm(ir v/w) 
indicated by vertical lines, the amplitudes diverge. For 
R c > R the amplitudes of the bouncing orbits quickly 
approach their asymptotic (zero-field) value. The inset 
shows this convergence in a wider range of B . 



for two-dimensional systems by Ozorio de Almeida 
and Hannay [28|, more explicit calculations have 
been performed, for example, by Kus et al. [29] 
and Sieber [|30[| . The main idea is always to re- 
place the saddle-point approximation by a better 
adapted uniform integration. An application to the 
disk billiard has not been attempted here and will 
be the subject of further studies. 

For the interpretation of the shell structure, let 
us first look at the weak-field regime (R c > R). 
The amplitudes for zero field given in Eq. ([!(]) are 
proportional to v -1 / 2 , favoring orbits with a small 
number of bounces v. The dependence of the am- 
plitudes on the magnetic field as shown in Fig. 18 
indicates that in the region where the "-" orbits dif- 
fer significantly from the "+" orbits, the latter are 
negligible. These effects^ together strongly favor 
the (2, 1) and the (3, 1) _ orbits. They end up with 
comparable amplitudes. From this picture we ex- 
pect as the dominating feature of the level density 
a pronounced beating pattern from the interference 
of the diameter and the triangular orbit. This beat- 
ing pattern is well known for the zero-field case. In 
three-dimensional metal clusters, it is usually re- 
ferred to as supershell oscillations |3j|J^] Our de- 
scription suggests that this beating will survive in 
homogeneous magnetic fields up to a strength of 



13 The G dependence of F(G) also supports slightly this 
effect. 

14 In the 3D spherical cavity, the beat comes from the inter- 
ference of the triangle and the square orbits (see Ref. |15|] ). 



B = kR. This is indeed observed, as seen in Fig. 16 . 
The thick lines in the frames (la) and (lb) corre- 
spond to a function^ 



sin(fcG(i )2 )) + sin(A:G( li 3)- 

. f AG\ . (,G^ 

sin k sin k — 

2 I 2 



(17) 



It predicts correctly the structure of the level den- 
sity in this regime. 

Approaching the field strength where R c = R, 
all orbits change G sharply to 2tt, so that they in- 
terfere coherently, forming the whispering gallery 
mode. We therefore expect that the beating behav- 
ior will disappear, leaving just a simple oscillation 
with the common frequency. In Fig. |l6] this sudden 
stop of the beat at R c = R can be clearly seen. The 
solid line in frame 2 shows that the frequency of the 
remaining single oscillation is predicted correctly. 

For R c < R, the influence of the cyclotron or- 
bits increases with stronger fields. The large am- 
plitudes of the bouncing orbits near the bifurcation 
points is, as we have already pointed out, unphys- 
ical and should be removed by a rigorous treat- 
ment of the orbit bifurcations. For strong fields, 
only cyclotron orbits and bouncing orbits with a 
great number of bounces v exist. The amplitudes 
of the latter are proportional to v^ 1 ^ 2 , so that in 
very strong fields we expect that the cyclotron or- 
bits dominate the level density. The gray lines in 
frame 3 of Fig. |l^ show the corresponding oscillat- 
ing term ,[^1 which, indeed, reproduces the main fea- 
ture of the quantum- mechanical result (solid black) . 
The skipping orbits with greatest amplitudes are 
those close to their bifurcation points. As can be 
seen in Fig. |l7|, all those orbits have nearly the 
same value of G = w ■ ir 2 R. Their contributions 
should therefore interfere constructively, giving rise 
to small structures in the level density of this pe- 
riod. Such structures can indeed be observed in 



a higher-resolution spectrum, as shown in Fig. |15. 
The spacing of the small peaks between the Lan- 
dau levels is, indeed, consistent with our simple pic- 
tureQ 



15 The phases are, of course, adjusted. 

16 For a simpler comparison, the amplitude is chosen to rise 



quadratically, as indicated by Eq 

17 This holds for the spacing of levels that "belong" to the 
same Landau level, and as long as we still have skipping 
orbits and do not enter the grazing orbit regime, where the 
reflection phases change. 
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Altogether we could show that this simple semi- 
classical picture is able to explain the main features 
of the quite complicated behavior of the level den- 
sity for arbitrarily strong fields in terms of just 3 
classical periodic orbits. We have here interpreted 
the dependence of the level density on the energy, 
but a completely analogous approach for the depen- 
dence on the magnetic field is possible. 

6 Summary 

We have derived a trace formula for the oscillating 
part of the level density of a circular billiard in ho- 
mogeneous magnetic fields. We have used the gen- 
eral approach of Creagh and Littlejohn and com- 
pared our findings with the quantum-mechanical 
solution. In the weak-field domain, where the clas- 
sical cyclotron radius R c is larger than the disk ra- 
dius R, the agreement is excellent and even a full 
quantization, e.g., the resolution of the level density 
into individual energy levels, is possible. 
In stronger fields, the quality of the standard semi- 
classical approximation is not satisfactory, even for 
the gross-shell structure. We have identified bound- 
ary effects to be responsible for the major part of 
the deviations. To implement these effects in the 
semiclassical trace formula, we have replaced the 
(discrete) Maslov index by a (continuous) reflec- 
tion phase. The latter was calculated in a simple 
one-dimensional approximation. With this correc- 
tion, the semiclassical approximation to the exact 
quantum-mechanical level density is good for all 
field strengths and energies. For a correction of 
the remaining deviations it would be necessary to 
include diffractive orbits and the effects of the or- 
bit bifurcations. The orbit bifurcations at strong 
field strengths affect the shell structure only to a 
small extent, their influence on the full quantiza- 
tion is even smaller. The diffractive orbits do not 
influence the shell structure but only the full quan- 
tization. Both effects can therefore be neglected for 
the semiclassical description of the gross-shell struc- 
ture. The reflection phases, however, are a crucial 
correction for both the gross-shell structure and the 
full quantization. 

One advantage of the semiclassical description 
is its easy numerical evaluation. Much more at- 
tractive, however, is the simple, intuitive picture 
gained from it. Quantum mechanics readily gives 
information on individual levels or level statistics, 
which are hard to derive semiclassically. But the ex- 



perimentally important long-range correlations of 
levels, leading to shells and supershells, are very 
easy to explain semiclassically. For a qualitative 
description of the shell structure just one or two 
classical periodic orbits are sufficient. In strong 
fields the single oscillation of the cyclotron orbits 
dominates and the coherent superposition of the 
strongest skipping orbits gives rise to additional 
small structures with much smaller spacing. For 
field strengths with R C <R the skipping orbits form 
coherently the whispering gallery mode, which gives 
rise to a single oscillation of the level density. In 
weak fields, the interference between the diameter 
and the triangular orbit dominates the level density. 
A quantitative semiclassical description is already 
possible including between 10 and 20 orbits. Fu- 
ture studies will be aimed at a rigorous treatment 
of the orbit bifurcations and an implementation of 
diffractive effects. 
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A Geometrical quantities of the 
PO 

The geometrical lengths c, d, and s and the an- 
gles Q, 7, and ip sketched in Fig|] can be expressed 
in terms of the classical cyclotron radius R c , the 
disk radius R and the classification parameter (3^ = 
(v, 11))^ as follows: 

G = — 7T , 

V 

■ ( R ■ n\ 

7 = arcsml— sinH I , 

f 7 - G + vr/2 for , R c > R) 
if = I -7 + G + vr/2 for (/3 + ,i? c < R) 
{ 7 + G - it/2 for (/T) , 

c = R cos <p , 

s = \]R C 2 -R 2 sin 2 G , 
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\s — i?cos 6 
s + R cos 8 



for ^ 
for p 



(18) 



B Evaluating the PO sum 

Semiclassical trace formulas are asymptotic series 
with non trivial convergence properties, so that they 
cannot be summed up straightforwardly. Frequently 
the Gaussian smoothing technique is used, which 
approximates the level density folded with a Gaus- 
sian by the trace formula where the amplitudes are 
damped by an additional (Gaussian) factor. This 
approach is limited to Gaussian line shapes and 
to smoothing of the level density in k. In this 
Appendix we introduce a more general approach, 
which can deal with arbitrary line shapes and smooth- 
ing variables. We will also state explicitly the con- 
ditions for the approximation to be valid. 



The general form of a trace formula is given by 

(19) 



^ — _ . 3-p (E) . 77- 

5g = J2ME)e 1 ^-^ 
r 



where T is a one-dimensional classification of the 
classical periodic orbits. If there is a generalized en- 
ergy e(E), and functions G(T,E) and <r(G), which 
fulfill 



- a r - -eG- o(G) 
we can rewrite the trace formula as 

5g = ^2A 2 (e,G)e ieG . 
G 



(20) 



(21) 



Rescaling G we can always obtain G G IN; the 
rescaling factors should be included in A 2 (e, G). Let 
us first assume that A 2 factorizes in terms that only 
depend on the generalized energy e and the classi- 
fication variable G: 

A 2 (e,G) = A G (G) A e (e) . (22) 

Approximating Eq. (|2l]) by an integral 

5g « A e (e) [ A G (G)e ieG dG . (23) 

JG 

gives (apart from normalization constants) the os- 
cillating part of the level density 5g as the Fourier 
transform of Ag(G): 

6g(e) « Vto A e {e) F[Ag(G)] . (24) 



For an arbitrary window function F(G) we get, us- 
ing the well-known folding theorem, 

/ F(G)A 2 (e,G)e ieG dG^5g(e)*f(e) . (25) 
J G 

Here /(e) denotes the Fourier transform of F(G) 
and "*" stands for the convolution integral. There- 
fore we have 



grCg) 



-MTp- 



Sg(e)* /(e) 



5g F :=Y / F(G)A r (E) e 

(26) 

where 5g F denotes the trace formula with damped 
amplitudes. This relation shows that folding the 
semiclassical level density with a smoothing func- 
tion /(e) is equivalent to a multiplication of the 
amplitudes with a window function F(G). Unfor- 
tunately the restrictions of Eqs. ( p0| ) and ( [22] ) are 
quite severe and often prevent the application of 



Eq. (26). With two additional approximations we 
can relax these restrictions. In the general case 
Eq. (p2[ ) is violated and we may just separate out a 
common dependence of the amplitudes on e: 



A 2 (e,G)=A G (e,G) A e (e) 



(27) 



In this case Eq. (|26D is still a good approximation 
if Ac(e, G) is sufficiently slowly varying in e. If we 
denote the characteristic width of /(e) with 7, this 
means that Aa(e, G) has to be nearly constant over 
a region 7 in e. If, on the other hand, there are no 
functions e(E) and G(E,F) that fulfill Eq. @, a 
local expansion of the action S in powers of e can 
be used: 



S 

h 



S(e ) 

h 



+ G(e ) (e-e )+O(e-e ) i 



(28) 



If this approximation is good in a region in e wider 
than the typical width 7 of the smoothing func- 
tion, Eq. (f26"l) still holds. In the general case G is 
therefore given by the first derivative of the classical 
action with respect to e: 



G(E) 



IdS 

h de 



(29) 



With e = E, fiG is the period T of the orbit, so 
that we refer to TiG as the quasiperiod. For systems 
with constant absolute velocity along the orbit (this 
holds especially for billiards), we get for the choice 
e = k 



dS 



dS dE _ kfi 2 
de dE dk m 
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where L is the geometrical orbit length. Putting 
all approximations together, we have shown that 
damping the amplitudes in the trace formula with 
a window function depending on G gives an approx- 
imation for the level density folded with the Fourier 
transform of the window function used: 



are Eqs. ( |3l| ) and (^2|), which give the limits of va- 



$g F ~ /(e) * 5g 



(30) 



This is the main result of this Appendix. The ap- 
proximation holds if in a region wider than the typ- 
ical width 7 of the smoothing function the condi- 
tions 

S « 5(e ) + G(e ) (e - e ) (31) 



and 



Azie, G) const 



(32) 



are fulfilled. These conditions depend mainly on 
the behavior of the actions and amplitudes. In or- 
der to match them, a well-adapted choice of the 
generalized energy is essential. Note that for nar- 
row smoothing functions (small 7), the conditions 
are less restrictive. Therefore for a full quantiza- 
tion the use of Eq. ( |3Q|) is often justified, whereas 
for the calculation of the gross-shell structure the 

put tight limits on 
which 



conditions Eqs. (31) and (| 
the use of the amplitude damping ansatz 
might seem counter-intuitive at first sight. 

We now illustrate the result with a simple exam- 
ple. Pure billiard systems are those where the the 
action along the orbits scales with the wave num- 
ber: S = hk ■ L, and L, the geometric orbit length, 
is independent of the energy. Setting 



e(E) = k 



<2mE 



and G(T) = L 



(33) 



Eq. (§ID is fulfilled trivially. If Eq. @ is also 
matched, then the use of a window function F de- 
pending on the orbit length L is equivalent to a fold- 
ing of the level density in k. Using a Gaussian win- 
dow function we get a Gaussian smoothing of the 
level density in k space. This is the technique fre- 
quently applied when evaluating trace formulas for 
billiard systems. Equation ( |30[) is somewhat more 
general, as it is not restricted to billiard systems 
nor to special window functions. It makes (at least 
in principle) the calculation of arbitrary line shapes 
within the POT possible. It can also be used for an 
estimation of the effects of a (numerical) truncation 
of the trace formula, which can be thought of as a 
special window function. More important, however, 



lidity of the amplitude damping formula (3C). 



B.l Evaluation for the circular billiard 

We want to apply the considerations of the last sec- 
tion on the circular billiard. The natural choice for 
the generalized energy is k. Then the quasiperiod 
is the geometrical orbit length, given by 

G = vr J 2. -2 7 for W+ : R c <R) 
27 otherwise. 

Note that for R c > R (weak fields) G is indepen- 
dent of the direction of motion ±. 
For computing the trace formula we have to choose 
an appropriate window function. As we want to 
compare the semiclassical result with the exact quan- 
tum-mechanical one, we look for a window function 
that can be Fourier transformed analytically. The 
usual Gaussian is nonzero for all G and has to be 
truncated, being thus no longer analytically Fourier 
transformable. We used a triangular window in- 
stead, which matches all our demands. In order to 
make our results comparable with the usual Gaus- 
sian smoothing, we characterize the window func- 
tion with a parameter 7, which corresponds to the 
variance of a Gaussian exp[— l/2(A;/7) 2 ] with the 
same half-width. 



We still have to check if the conditions (31) 
and ( |32"|) hold. They depend on the behavior of the 
amplitudes that are plotted in Fig. At the bifur- 
cation points the orbit amplitudes diverge, so that 
Eq. (32) is violated. For the evaluation in the corre- 
sponding regions we have therefore used a numer- 
ical folding procedure and evaluated directly the 
right-hand-side of Eq. (^0[) . 

For the cyclotron orbits discussed in Sec. |4.3| we get 
G = n- 2ttR c and A = {2E Q )- 1 (l - R c /R) 2 , which 
is slowly varying in the whole energy range. For 
the cyclotron orbits, approximation ( pOj ) is there- 
fore justified for all E and B. 
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